Circulating metabolic profile in idiopathic pulmonary fibrosis: data from the IPF-PRO Registry

Background The circulating metabolome, reflecting underlying cellular processes and disease biology, has not been fully characterized in patients with idiopathic pulmonary fibrosis (IPF). We evaluated whether circulating levels of metabolites correlate with the presence of IPF, with the severity of IPF, or with the risk of clinically relevant outcomes among patients with IPF. Methods We analyzed enrollment plasma samples from 300 patients with IPF in the IPF-PRO Registry and 100 individuals without known lung disease using a set of targeted metabolomics and clinical analyte modules. Linear regression was used to compare metabolite and clinical analyte levels between patients with IPF and controls and to determine associations between metabolite levels and measures of disease severity in patients with IPF. Unadjusted and adjusted univariable Cox regression models were used to evaluate associations between circulating metabolites and the risk of mortality or disease progression among patients with IPF. Results Levels of 64 metabolites and 5 clinical analytes were significantly different between patients with IPF and controls. Among analytes with greatest differences were non-esterified fatty acids, multiple long-chain acylcarnitines, and select ceramides, levels of which were higher among patients with IPF versus controls. Levels of the branched-chain amino acids valine and leucine/isoleucine were inversely correlated with measures of disease severity. After adjusting for clinical factors known to influence outcomes, higher levels of the acylcarnitine C:16-OH/C:14-DC were associated with all-cause mortality, lower levels of the acylcarnitine C16:1-OH/C14:1DC were associated with all-cause mortality, respiratory death, and respiratory death or lung transplant, and higher levels of the sphingomyelin d43:2 were associated with the risk of respiratory death or lung transplantation. Conclusions IPF has a distinct circulating metabolic profile characterized by increased levels of non-esterified fatty acids, long-chain acylcarnitines, and ceramides, which may suggest a more catabolic environment that enhances lipid mobilization and metabolism. We identified select metabolites that were highly correlated with measures of disease severity or the risk of disease progression and that may be developed further as biomarkers. Trial registration ClinicalTrials.gov; No: NCT01915511; URL: www.clinicaltrials.gov. Supplementary Information The online version contains supplementary material available at 10.1186/s12931-023-02644-7.


Background
Idiopathic pulmonary fibrosis (IPF) is a progressive and ultimately fatal interstitial lung disease [1].IPF is characterized by progressive decline in lung function due to the excess deposition of extracellular matrix components, which ultimately leads to destruction of the lung architecture and respiratory failure [2].Epithelial cell injury with activation of fibroblasts and production of pro-fibrotic growth factors is considered pathogenic [3].Two antifibrotic drugs, nintedanib and pirfenidone, slow decline in lung function in patients with IPF [4,5], and have been licensed for its treatment, but the disease remains progressive.The rate of disease progression among patients with IPF is variable [6,7] and while some circulating biomarkers have shown promise as predictors of progression [8][9][10][11], for an individual patient, it remains challenging to predict the rate at which their disease will progress.
As metabolites represent the end products of cellular processes and are potential effectors of disease biology, the circulating metabolome may provide novel insights into disease activity and identify candidate biomarkers of disease progression.This idea is supported by studies suggesting that the pathogenesis of pulmonary fibrosis is associated with metabolic abnormalities including mitochondrial dysfunction; lipid dysregulation; elevations in lactic acid and lactate dehydrogenase 5; alterations in glycolysis, glutathione biosynthesis, adenosine triphosphate degradation, and ornithine aminotransferase pathways [12][13][14][15][16][17][18][19][20].The circulating metabolome in patients with IPF has not been extensively studied, but associations have been identified between select metabolites and measures of disease severity and the risk of progression [20,21].In this study, we quantified the circulating metabolome in patients with IPF and controls without known lung disease to identify metabolites present at different levels between the groups.We then related circulating metabolite levels to measures of disease severity in the patients with IPF and evaluated their utility for predicting the risk of clinically relevant outcomes among these patients.

Study population
The IPF cohort was drawn from the Idiopathic Pulmonary Fibrosis PRospective Outcomes (IPF-PRO) Registry, a multicenter observational registry of patients with IPF that was diagnosed or confirmed at the enrolling center in the past 6 months [22].The cohort comprised 300 patients, enrolled between June 2014 and February 2017, who had a blood sample and data on critical clinical variables at enrollment, including diagnostic criteria for IPF (definite, probable, possible) according to the 2011 international guidelines [23].
The control cohort was drawn from the Measurement to Understand the Reclassification of Disease of Cabarrus/Kannapolis (MURDOCK) Study, a registry of adult residents in North Carolina in which self-reported health information and biological samples are collected [24].To ensure the control cohort had a similar age, race, and ethnicity distribution to the IPF cohort, controls were White, non-Hispanic and aged 60 to 80 years.Participants were excluded if they had self-reported respiratory disease, cancer, or autoimmune disease, were active smokers, had active second-hand tobacco exposure, or used respiratory-targeted medications or immunomodulators.Random sampling with stratification by sex and smoking status (ever versus never) was used to select 100 controls with a similar distribution of these characteristics to the IPF cohort.
The IPF-PRO Registry was approved by the Duke University Institutional Review Board (Pro00046131).The IPF-PRO Registry protocol was also approved by the relevant Institutional Review Boards and/or local Independent Ethics Committees prior to enrollment at each site (listed in the Acknowledgments).All patients provided written informed consent.The IPF-PRO Registry was registered with ClinicalTrials.gov(Identifier: NCT01915511).The MURDOCK Study Community Registry and Biorepository was approved by the Duke University Health Institutional Review Board (Pro00011196) and all participants provided written informed consent.The MURDOCK Study Community Registry and Biorepository was registered with Clinical-Trials.gov(Identifier: NCT01708408).

Metabolite quantification
A set of metabolites comprising 15 amino acids, 45 acylcarnitines, 21 ceramides, 34 sphingomyelins, 3 branched chain keto-acids, and 3-hydroxyisobutyrate was quantified in plasma samples taken at enrollment using flowinjection tandem mass spectrometry (MS/MS) [25,26] or liquid chromatography(LC)-MS/MS as previously described [27][28][29][30].Quantitative measurements were achieved by adding known quantities of stable isotopelabeled internal standards to the biological samples.The assays were run in a 96-well-plate format, with a calibration curve and a set of two quality-control samples at the beginning and end of each plate.An ion ratio of the analyte respective to the internal standard was computed from centroided spectra, and ion ratios converted to concentrations using the calibration curve slope; this allowed values to be analyzed even if below the limit of quantification.Values of 0 were imputed as half the lowest detectable value for a given analyte.
To provide context for interpretation of the metabolite data, several common clinical analytes were measured, using published methods [28,29], in each sample: total, high-density lipoprotein (HDL) and low-density lipoprotein (LDL) cholesterol, glucose, ketones, lactate, non-esterified fatty acids (NEFA), triglycerides, glycerol, and 3-hydroxybutyrate.There was insufficient sample volume to permit measurement of 3-hydroxybutyrate in 9.2%, HDL cholesterol in 4.0%, LDL cholesterol in 1.0% and ketones in 0.2% of samples.For analyses of clinical analytes, only cases and controls with complete data for all analytes were included.Metabolite and clinical analyte concentrations were log 2 transformed prior to analysis.

Statistical analyses
Descriptive statistics were used to analyze patient characteristics and the level of each metabolite or clinical analyte in patients with IPF and controls.Linear regression was used to assess whether metabolite or clinical analyte concentrations differed by IPF or control status.P values were corrected for multiple comparisons using the Benjamini-Hochberg procedure to control the false discovery rate (FDR) at 5%.Comparisons were considered statistically significant if the FDR-corrected P < 0.05 and clinically significant if the absolute fold-change ≥ 30% (i.e.absolute log 2 fold-change ≥ 0.38).
Among patients with IPF, linear regression models were used to determine the associations between circulating metabolites and three clinical measures of disease severity at enrollment: forced vital capacity (FVC) % predicted, diffusion capacity of the lung for carbon monoxide (DLco) % predicted, and composite physiologic index (CPI) [31].Each measure was analyzed as a continuous variable.Comparisons were considered statistically significant if the FDR-corrected P < 0.05 and clinically significant if ≥ 5-unit difference in the disease severity measure per unit change in the log 2 uM concentration of the metabolite (i.e. if a doubling of the metabolite concentration was associated with a ≥ 5-point difference in the disease severity measure).
To identify candidate biomarkers of clinically relevant outcomes among patients with IPF, univariable associations between levels of metabolites or clinical analytes at enrollment and outcomes were determined using Cox proportional hazards regression analyses.We present results for all-cause death; respiratory death; composite of all-cause death or lung transplantation; composite of respiratory death or lung transplantation; composite of decline in FVC ≥ 10% predicted, death, or lung transplantation; composite of decline in DLco ≥ 15% predicted, death, or lung transplantation.Kaplan-Meier plots were used to describe cumulative event probabilities.Analyses were unadjusted and adjusted for clinical factors known to influence outcomes i.e., sex, age, FVC % predicted, DLco % predicted, supplemental oxygen use (all assessed at enrollment).Comparisons were considered statistically significant if the FDR-corrected P < 0.05 and clinically significant if the hazard ratio was < 0.67 or > 1.5.The linearity and proportional hazards assumptions were assessed.For metabolites for which the linearity assumption failed, piece-wise linear splines with 1 or 2 knots were used to characterize the non-linearity and hazard ratios with 95% confidence intervals for each segment are presented to describe the relationship between the metabolite and outcome.For metabolites for which the proportional hazards assumption failed, an interaction term with time was included in the model as a time-dependent covariate and hazard ratios with 95% confidence intervals at 6, 12, and 24 months are presented to describe how the association changed during follow-up.

Associations between metabolites and IPF
The concentrations of metabolites and clinical analytes in the IPF and control cohorts are given in Tables S1 and S2 in Additional file 1. Linear regression analyses identified 64 metabolites and 5 clinical analytes with significantly different levels between patients with IPF and controls.Sixteen of these had an absolute fold-change ≥ 30% (log 2 fold-change ≥ 0.38) between groups and four had an absolute fold-change ≥ 50% (log 2 fold-change ≥ 0.58) between groups (Fig. 1; Tables S3 and S4 in Additional file 1).Levels of these metabolites/clinical analytes were not significantly different between patients with IPF who were versus were not using anti-fibrotic therapy at enrollment (Table S5 in Additional file 1).Among the clinical analytes and metabolites with the greatest difference between groups were NEFA, the amino acid ornithine, multiple long chain acylcarnitines and select ceramides, with higher levels seen in the patients with IPF relative to the controls (Fig. 1).

Circulating metabolites and outcomes
Outcomes were assessed over a median (25th, 75th percentile) follow-up period of 39.9 (20.1, 53.3) months.The cumulative event probability of each outcome is shown in Fig. 2. In unadjusted analyses, select amino acids, acylcarnitines, ceramides and sphingomyelins were significantly associated with outcomes (Fig. 3).After adjusting for clinical factors, higher levels of the acylcarnitine C:16-OH/C:14-DC remained significantly associated with all-cause mortality, lower levels of the acylcarnitine C16:1-OH/C14:1DC remained significantly associated with all-cause mortality, respiratory death, and respiratory death or lung transplant, and higher levels of the sphingomyelin d43:2 were associated with the risk of respiratory death or lung transplantation (Fig. E1 in Additional file 1).

Discussion
We found that the peripheral blood metabolome of patients with IPF was distinct from that of individuals without known lung disease of similar age and sex distribution.In addition, select metabolites were highly correlated with measures of disease severity and with the risk of clinically relevant outcomes, even after accounting for clinical factors known to impact such outcomes.In our study, the patients with IPF had significantly higher circulating levels of several lipid species, including multiple long-chain acylcarnitines and ceramides and one sphingomyelin species, compared with the controls.These findings are consistent with data showing dysregulated fatty acid or lipid metabolism in mouse models of IPF [16,32], elevated free fatty acids or dysregulation of fatty acid metabolism-related genes in the lungs of patients with IPF [16,33], and elevated circulating levels of fatty acids in patients with IPF [19,21,34].The pattern Fig. 1 Associations between circulating metabolites and clinical analytes in patients with IPF vs. controls of higher levels of NEFA along with higher levels of acylcarnitines and select ceramides in patients with IPF in our study points to enhanced availability of fatty acids, possibly due to lipolysis related to wasting.An increase in lipolysis is supported by the observation that mobilization of fatty acids from adipose triacylglycerol stores is accompanied by release of the glycerol backbone upon which the fatty acids were esterified, and glycerol levels were higher in the plasma of patients with IPF compared with controls.Higher NEFA levels may also be a marker of impaired action of insulin to suppress lipolysis (insulin resistance) [35], although this seems less likely because the patients with IPF had lower circulating glucose levels and a smaller proportion had diabetes than the controls.As anti-fibrotic treatments are associated with gastrointestinal side effects [4,5], we considered whether the derangements in circulating lipids might be related to weight loss and/or malabsorption, but among the lipid metabolites with higher levels in patients with IPF, we did not find differences in levels between patients with IPF who were and were not using anti-fibrotic therapy.Further evaluation of the relationships between weight change and circulating lipids in patients with IPF is required.
A study in patients with chronic obstructive pulmonary disease demonstrated inverse relationships between individual sphingolipids and the severity of emphysema, and measurements of select sphingolipids improved prediction of severe exacerbations beyond predictions based on demographic and clinical covariates [36].In studies in patients with IPF, negative correlations were observed between levels of select signaling lipids in peripheral blood and lung function [34], and higher levels of certain triglycerides and phosphatidylcholines were found in the plasma of patients with more rapid progression [21].A risk model based on expression of select fatty acid genes in bronchoalveolar lavage fluid predicted outcomes among patients with IPF [33].In our study, higher levels of ceramide C18:1/C16:0 was associated with more severe disease.Interestingly, this ceramide has been associated with inhibition of mitochondrial fatty acid oxidation [37], which may support the acylcarnitine profile described in our study in the sense that limits of mitochondrial oxidative activity could lead to accumulation of partially oxidized fatty acyl CoAs and acylcarnitines [38].Of note, this particular ceramide has been linked to defects in mitochondrial fragmentation and a plethora of metabolic disease phenotypes [39], supporting the concept of IPF as a systemic metabolic disorder.
We found evidence of amino acid dysregulation in patients with IPF.Specifically, patients with IPF had significantly higher circulating levels of ornithine compared with controls.Ornithine, a non-proteinogenic amino acid, is derived from the actions of arginase and serves as an important intermediate in the urea cycle and as a precursor to other metabolites, including citrulline, proline and various polyamines important for wound healing and cell proliferation [40].Of note, proline and its derivative hydroxyproline comprise approximately 23% of the amino acid content of collagen [41], suggesting that higher levels of ornithine in the blood may be reflective of enhanced collagen production in the lung.A previous study reported elevations in 4-hydroxyproline levels and of the polyamines putrescine and spermidine in the lungs of patients with IPF [42], supporting the concept that ornithine metabolism is dysregulated in this disease.Our study also indicated that higher levels of certain amino acids were correlated with less severe disease: higher levels of the branched chain amino acids valine, leucine and isoleucine were associated with higher DLco % predicted and lower CPI.In addition to being key substrates for energy metabolism and protein synthesis, these amino acids play key signaling roles in regulating growth and energy pathways.Indeed, they are major activators of the mammalian target of rapamycin (mTOR), a pro-growth enzyme that is hyperactivated in fibroblasts and epithelial cells in IPF [43][44][45].Whether higher circulating levels of these amino acids correlate with mTOR activity in the lung is a concept worthy of further investigation.
Our study has several strengths, including the multicenter nature of the registries from which the IPF and control cohorts were derived and the comparison of patients with IPF to control participants with similar distributions of age, sex, and smoking history.There are also some inherent limitations in our approach.Our cohort is a US-based population of predominantly White patients and its generalizability to other populations of patients with IPF is uncertain.Although we characterized a broad array of metabolites, our approach was targeted rather than discovery-based, so metabolites of potential importance may have been missed.The analyses of clinical analytes did not include all the participants in the study due to insufficient sample volumes.

Conclusion
The results of this study suggest that circulating metabolites may hold value as diagnostic, disease activity, or prognostic biomarkers for IPF and may provide insights into the pathobiology of this disease.These data provide support for the development of biomarker-inclusive algorithms that enable risk stratification for patients with IPF.Multicenter cohort studies such as the IPF-PRO Registry, with well-annotated serial clinical data and biological samples, are of great value for identification and validation of biomarkers to improve the management of IPF.
A podcast describing the key data presented in this manuscript is available at: https://www.usscicomms.com/respiratory/Todd/IPF-PRO-BIOmetabolomics.

Fig. 2 Fig. 3
Fig. 2 Cumulative probability of an event of outcomes over follow-up

Table 1
Clinical characteristics of the IPF and control cohorts at enrollment

Table 2
Association of circulating metabolites or clinical analytes with disease severity measures at enrollment in patients with IPF in unadjusted analyses.Analytes meeting statistical significance (FDR-corrected P < 0.05) or clinical significance (≥ 5-unit difference in disease severity measure per unit change in log 2 uM of the metabolite concentration) thresholds are shown Definition of abbreviations: FDR false discovery rate, DLco diffusion capacity of the lung for carbon monoxide, FVC forced vital capacity, CPI composite physiologic index.*Observed unit difference in the disease severity measure per unit change in log 2 uM of metabolite concentration (e.g.doubling of metabolite concentration)